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M. Peimbert 1 and A. Peimbert 1 , 
RESUMEN 

Presentamos una revision sobre la relevancia del artfculo de Peimbert y Costero (1969), en la determinacion de 
la composition qui'mica de regiones H II. Analizamos las evidencias observacionales en favor de la presencia de 
variaciones de temperatura en nebulosas gaseosas. Hacemos una breve mention a los metodos para tomar en 
cuenta a los iones no observados en la determinacion de las abundancias de los elementos. 

ABSTRACT 

We present a review about the relevance of the paper by Peimbert and Costero (1969), on the chemical 
abundance determinations of H II regions. We analize the observational evidence in favor of the presence of 
temperature variations inside gaseous nebulae. We make a brief mention of the methods used to estimate the 
contribution of the unobserved ions to the total chemical abundances. 
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1. INTRODUCTION 

The paper by M. Peimbert and R. Costero (1969, 
hereafter PC) presents the determination of the 
chemical composition of three galactic HII regions, 
the Orion nebula, M8 and M17. To derive the chemi- 
cal abundances PC took into account the presence of 
temperature variations inside gaseous nebulae based 
on the definitions of the average temperature, To, 
and the mean square temperature variation, t 2 , in- 
troduced by Peimbert (1967), that are given by 



T (V,i) 



J T e N e N,dV 
jN.N.dV ' 



and 



t 2 = 



J(T e - T ) 2 iV e iW 
T 2 / N e N t dV ' 



(1) 



(2) 



respectively, where T e , N e , and Ni are the electron 
temperature, the electron density, and the ion den- 
sity of the volume element, and V is the observed 
volume. 

PC includes simple equations to determine the 
chemical abundances in the presence of temperature 
variations, additional and more detailed equations to 
determine the chemical abundances have been pre- 
sented by Peimbert et al. (2002), Ruiz et al. (2003) 
and Peimbert et al. (2004). 
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Accurate abundances are needed to test models 
of: gaseous nebulae, stellar evolution, galactic chemi- 
cal evolution, and the observable universe as a whole. 
The presence of temperature variations produces bi- 
ases in the abundance determinations derived from 
collisionally excited lines that are typically in the 0.1 
to 0.5 dex range. PC provided significant advances 
in the procedures to determine abundances of H II 
regions and planetary nebulae. 

We will mention the observational evidence in fa- 
vor of large temperature variations based on seven 
independent methods, and some of the most relevant 
results that abundances, derived including the effect 
of temperature variations, have had in different ar- 
eas of astronomical research. We will define those 
gaseous nebulae with t 2 values < 0.01 as small tem- 
perature variation objects, and those with t 2 values 
> 0.01 as large temperature variation objects. 

Most papers dedicated to the determination of 
abundances in gaseous nebulae assume that temper- 
ature variations are small, and consequently that the 
abundances can be derived under the assumption of 
constant temperature. In section 2 we present evi- 
dence in favor of large temperature variations inside 
most gaseous nebulae. In section 3 we present abun- 
dance determinations that agree with the presence 
of large temperature variations based on: the com- 
parison of stellar abundances with gaseous nebular 
abundances, the comparison of solar and H II region 
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abundances within the framework of chemical abun- 
dance models of the Galaxy, and the comparison of 
the primordial helium abundance derived from H II 
regions compared with that derived from the Stan- 
dard Big Bang Nucleosynthesis. In section 4 we dis- 
cuss briefly the ionization correction factor method 
to derive total chemical abundances presented in PC. 

2. SIMULTANEOUS DETERMINATIONS OF T 
AND t 2 

To determine To and t we need to combine elec- 
tron temperatures, T e , based on two different meth- 
ods: one that weighs preferentially the high tem- 
perature regions and one that weighs preferentially 
the low temperature regions. Photoionization mod- 
els of chemically homogeneous HII regions predict 
t 2 values typically in the 0.003 to 0.01 range, while 
observations yield t 2 values typically in the 0.02 to 
0.06 range. 

These results imply that there are additional 
heating and cooling processes in H II regions that are 
not yet included in photoionization models. Many 
possible causes of temperature variations have been 
discussed in the literature some of them are: X- 
ray heating, cosmic-ray heating, density variations, 
deposition of mechanical energy, deposition of mag- 
netic energy, presence of shadowed regions, chemical 
inhomogeneities, dust heating, and transient effects 
due to changes in the ionizing flux. The source of 
large temperature variations in gaseous nebulae is 
still an open problem, and the relative importance 
of the effects causing the variations might be differ- 
ent for each nebula. 

In what follows we will mention some relevant 
papers that, based on seven independent methods, 
indicate the presence of large temperature variations 
in gaseous nebulae. 

2.1. Temperatures derived from the Balmer 
continuum and Balmer line intensities together 
with temperatures derived from collisionally 
excited line intensities 

The large temperature variations idea was based 
mainly on three results: a) the smaller T(Bac) val- 
ues than those derrived from T([0 III]) and T([N II]) 
(PC and Peimbert 1967) and the t 2 values around 
0.04 derived from the photoionization models by 
Hjellming (1966), that included only N, O and Ne 
as potential coolants of the model H II regions. In 
the seventies and eighties most observers obtained 
t 2 values from T(Bac) and T(0 III) of about 0.02 
± 0.04 that were consistent with t 2 = 0.00. Moreover 



photoionization models that included more chemi- 
cal elements than those considered by Hjellming pre- 
dicted t 2 values of about 0.005. Probably these two 
results together with the simpler equations to deter- 
mine chemical abundances led most workers to adopt 
t 2 = 0.00 and the temperature provided by T(0 III) 
to determine the abundances of gaseous nebulae. 

The accuracy of the temperature determinations 
improved and Liu and Danziger (1993) by combining 
T(Bac) with T(0 III) determinations found large t 2 
values for 14 planetary nebulae with an average value 
around 0.03 and errors considerably smaller than the 
t 2 values. Similarly Garcia-Rojas, Esteban and col- 
laborators (Garcia-Rojas et al. 2004, 2005, 2006, 
2007, Esteban et al. 2004) by combining T(Bac) with 
T(0 III) and T(0 II) determinations for eight galac- 
tic HII regions found t 2 values in the 0.01 to 0.056 
range with an average value of 0.034, again with er- 
rors smaller than the t 2 values. Moreover Gonzalez 
Delgado et al. (1994) studied the giant extragalactic 
HII region NGC 2363 and based on measurements 
of the Paschen discontinuity found t 2 values of 0.064 
and 0.098 for knots A and B respectively. 

2.2. Abundances derived from collisionally excited 

C III lines together with abundances derived 
from C II recombination lines 

The N(C+ + /N(U+) values derived for HII re- 
gions and planetary nebulae based on the C II A 
4267 to H/3 intensity ratio are, in general, higher than 
those derived from the C III AA 1906 + 1909 to H/3 
intensity ratio; in some cases the difference reaches 
a factor of ten. General discussions of this prob- 
lem have been given in the literature (e. g. Torres- 
Peimbert, Peimbert & Daltabuit 1980; Kaler 1986; 
Rola & Stasinska 1994, Peimbert, Torres-Peimbert, 
& Luridiana 1995, Liu 2006, Peimbert & Peimbert 
2006, and references therein). Several ideas have 
been advanced to explain the discrepancy: a) errors 
in the atomic data, b) errors in the observations and 
c) the presence of spatial temperature variations. 

2.3. Abundances derived from collisionally excited 
[O III] lines together with abundances derived 

from O II recombination lines 

Peimbert, Storey, & Torres-Peimbert (1993) were 
the first to determine O/H values for gaseous nebu- 
lae based on the recombination coefficients for O II 
lines computed by Storey (1994). The temperature 
dependence of the O II lines is relatively weak and 
very similar to that of the HI lines, therefore the 
++ /H + ratios are independent of the electron tem- 
perature. Alternatively the ++ /H+ ratios derived 
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Method 



TABLE 1 

t 2 VALUES DERIVED FROM DIFFERENT METHODS 

Orion Nebula 30 Doradus NGC 5315 NGC 6543 



T(Bac) and T e ([0 II], [O III]) 
T(Hc II) and T e ([0 II], [O III]) 
V(C++)RL and V(C++)CEL 
V(0+)RL and V(0+)CEL 
V(0++)RL and V(0++)CEL 
V(Ne++)RL and V(Nc++)CEL 
Average 



0.018±0.018 0.022±0.007 0.039±0.022 0.028±0.009 

0.022±0.002 ••• 0.060±0.007 0.035±0.014 

0.039±0.011 0.056±0.040 0.063±0.035 0.036±0.010 

0.052±0.029 0.075±0.040 

0.020±0.002 0.038±0.005 0.048±0.004 0.024±0.008 

0.032±0.014 ••• 0.068±0.020 0.022±0.010 

0.022±0.002 0.033±0.005 0.051±0.004 0.028±0.005 



from collisionally excited lines do depend strongly on 
the average temperature, To, and the mean temper- 
ature square, t 2 (e. g.: Peimbert 1967, PC, Ruiz et 
al. 2003, Peimbert et al. 2004). 

Garcfa-Rojas, Esteban and collaborators 
(Garcfa-Rojas et al. 2004, 2005, 2006, 2007, Este- 
ban et al. 2004) observed 8 galactic H II regions 
and by combining the [O III] lines with the O II 
lines find t 2 values in the 0.020 to 0.046 range 
with an average value of 0.038. Similarly, Esteban 
et al. (2002), Peimbert (2003), Bresolin (2007), 
and Esteban et al. (2009) for 11 extragalactic HII 
regions find t 2 values in the 0.027 to 0.124 range 
with an average value of 0.048. 

2.4. He I recombination lines 

From a large set of He I lines it is possible to de- 
termine T(He I). This temperature combined with 
T([0 III]) and T([0 II]) yields T and t 2 (Peim- 
bert, Peimbert, & Ruiz 2000, Peimbert, Peimbert, 
& Luridiana 2002). This method has been applied 
by Garcfa-Rojas, Esteban and collaborators (Garcfa- 
Rojas et al. 2004, 2005, 2006, 2007, Esteban et al. 
2004) to 7 galactic H II regions yielding t 2 values in 
the 0.017 to 0.046 range, with an average value of 
0.027. 

This method has also been applied by Peimbert, 
Luridiana, and Peimbert (2007a) to five metal poor 
extragalactic HII regions obtaining an average t 2 
value of 0.026; and by Esteban et al. (2009) to 14 
giant extragalactic H II regions obtaining t 2 values 
in the 0.022 to 0.125 range with an average value of 
0.058. 

2.5. Comparison of the t 2 values derived from 
different methods 

In Table 1 we present t 2 values determined from 
six different methods for four well observed objects, 
two H II regions and two planetary nebulae (Esteban 
et al. 2004, Peimbert 2003, Peimbert et al. 2004, 



Georgiev et al. 2008). For each object the different 
methods yield t 2 values that are in very good agre- 
ment. Moreover the adopted average values show er- 
rors that are from five to twelve times smaller than 
the t 2 determinations. These two results imply that 
the temperature variations are real and very large. 

To explain the presence of large spatial temper- 
ature fluctuations, considerably higher than those 
predicted by chemically homogeneous photoionized 
models, several possibilities have been proposed. Re- 
views on this problem have been presented by many 
authors sec for example: Torres-Peimbert, Peim- 
bert, & Peha (1990), Peimbert, Sarmicnto, & Fierro 
(1991), Esteban (2002), Torres-Peimbert and Pe- 
imbert (2003), Liu (2006), Peimbert, & Peimbert 
(2006). 

2.6. High spatial resolution map of the columnar 
temperature in the Orion nebula 

O'Dell, Peimbert, and Peimbert (2003), based on 
observations with the Hubble Space Telescope, pre- 
sented a high spatial resolution map of the columnar 
electron temperature T c of a region to the south west 
of the Trapezium in the Orion Nebula . From their 
T c (0 ++ ) map of the Orion Nebula, that includes 1.5 
x 10 6 independent temperature determinations, they 
found that the observed mean square temperature 
variation in the plane of the sky, t 2 a (0 ++ ), amounts 
to 0.008. 

Based on their t 2 a (0 ++ ) value, together with 
geometrical considerations and other observations in 
the literature, they estimated that t 2 (0 ++ )= 0.021. 
Note that the total t 2 (0++) is larger than t 2 a (0++) 
because in addition to the variations across the plane 
of the sky it includes the temperature variations 
along the line of sight. From their t 2 (0 ++ ) value 
and comparisons between the temperatures in the 
low- and high-ionization zones, the + and ++ 
zones, they found that t 2 (H n)= 0.028 ± 0.006. 
Their derived t 2 (H n) value is 7 times higher than 
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those obtained from homogeneous one-dimensional 
photoionization models of the Orion Nebula (e. g. 
Baldwin et al. 1991, Peimbert et al. 1993). 

2.7. Comparison of the different ways to calibrate 
Pagel's method 

The most popular metallicity indicator to de- 
termine the O /H ratio in extragalactic H II re- 
gions was introduced by Pagel et al. (1979, see 
also Edmunds & Pagel 1984) and is indistinctly 
known as Pagel's, or i?23> or O23 indicator, where 
023 = 7([Oll]A3727 + [OIII]AA4959,5007)/7(H/?). 
The O23 indicator has been calibrated with the 
O/H values based on three different methods: a) 
by using photoionization models, that we will call 
PIM calibrations or PIM method, b) by using ob- 
servational determinations of the O/H abundances 
based on the electron temperature derived from 
the 7(4363)/7(5007) [OIII] ratio together with the 
7(3727)/7(H/3) and the 7(5007) /7(H/3) line ratios, 
the so called T(4363) method, and c) by using ob- 
servational determinations of the O/H abundances 
based on the intensity ratio of O II recombination 
lines to H I recombination lines that has been called 
the Olliu, method by Peimbert et al. (2007b). 

These three methods have been discussed by sev- 
eral authors, (e. g. McGaugh 1991, Pilyugin & 
Thuan 2005, Peimbert et al. 2007b, Kewley & Elli- 
son 2009, and references therein). These three meth- 
ods produce very different O/H calibrations, in ex- 
treme cases the differences in the inferred O /H abun- 
dances among these calibrations reach values of 0.7 
dcx. 

Peimbert et al. (2007b) reach the following con- 
clusions: a) the O Hrl method supports the sugges- 
tion that the controversy produced by the relatively 
high O/H values predicted by the PIM calibrations 
and the relatively low O/H values predicted by the 
T(4363) calibrations are mainly due to temperature 
variations; b) the best way to calibrate the O23 in- 
dicator is to use the Hrl method to obtain the 
O/H values because it is independent of the tem- 
perature structure of the H II regions; c) the use of 
T(4363) values to derive O/H, under the assump- 
tion of constant temperature, provides a lower limit 
to the O/H abundance ratios; d) since the nebular 
lines are less sensitive to To and t 2 than the auroral 
lines, the model calibrations that adjust the nebu- 
lar lines are closer to the O IIrl calibration than 
those derived using the observed T(4363) values; 
and e) the presence of temperature variations af- 
fects strongly the T(4363) method, weakly the PIM 
method, and leave the O IIrl method unaffected, or 



TABLE 2 

STELLAR AND NEBULAR ABUNDANCES FOR 
NGC 6543 A 



Element 


Stellar 


Nebular(RC) Ncbular(CL) 


He 


11.00±0.04 


11.05±0.01 




C 


9.03±0.10 


8.90±0.10 


8.40±0.10 


N 


8.36±0.10 


8.83±0.20 


8.43±0.20 


O 


9.02±0.10 


9.19±0.12 


8.79±0.06 


s 


7.57±0.10 




7.08±0.06 


Ne 




8.67±0.10 


8.25±0.06 



a In units of 12 + Log N{X)/N(H). 



in other words the O II_rl method is independent of 
the temperature structure of the nebula. 



3. ADDITIONAL SUPPORT IN FAVOR OF 
LARGE TEMPERATURE VARIATIONS 

3.1. Comparison of the abundances of the nebula of 
NGC 6543 with those of its central star 

Table 2 presents the stellar abundances of the 
planetary nebula NGC 6543 together with the nebu- 
lar envelope abundances derived from recombination 
lines (RC) that do not depend on the temperature 
structure of the nebula, and nebular abundances de- 
rived from collisionally excited lines (CL) under the 
assumption that t 2 = 0.00 determined by Georgiev 
et al. (2008). 

The stellar abundances of He, C and O are in ex- 
cellent agreement with the nebular abundances de- 
rived from recombination lines. In particular the 
similar He/H values imply that there are no He rich 
inclusions present in the nebula. On the other hand 
the C, O and S abundances derived from nebular 
collisionally excited lines (under the assumption of 
t 2 = 0.00) are considerably smaller than the stellar 
abundances. Table 1 shows the t 2 values needed to 
reach agreement between the recombination and the 
collisional abundance determinations. These results 
imply that indeed large temperature variations are 
present NGC 6543 and that the nebula of this object 
is chemically homogeneous. 

In addition a t 2 value of about 0.028 permits to 
reach agreement between the nebular determination 
of S and the stellar one. Similarly a t 2 value of about 
0.028 also permits to reach agreement between the 
Ne abundances derived from recombination and col- 
lisionally excited lines. 
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3.2. Comparison of the oxygen abundances of the 
Orion nebula with those of B stars of the solar 

vicinity 

Przybilla, ct al. (2008) based on the study of 
B stars of the solar vicinity have found that 12 + 
log O/H = 8.76 ± 0.03, while Simon-Diaz (2009 pri- 
vate communication) finds that 12 + log O/H = 8.76 
± 0.04 for the Orion association B stars, these results 
are in excellent agreement with that derived by Es- 
teban et al. (2004) based on the O recombination 
lines, that amounts to 8.73 ± 0.03. The result by 
Esteban et al. (2004) includes a correction of 0.08 
dex due to the fraction of O tied up in dust grains 
estimated by Esteban et al. (1998). By adopting 
the O dust correction of 0.12 ± 0.03 dex estimated 
by Mesa-Delgado et al. (2009), the Orion nebula O 
abundance derived by Esteban et al. (2004) becomes 
8.77 ± 0.04, also in excellent agreement with those 
derived from the B stars. 

On the other hand, based on the T(4363) method 
with t 2 = 0.00 Dcharveng et al. (2000), Pilyu- 
gin, Ferrini, & Shkvarun (2003), and Esteban et al. 
(2004) obtain for the Orion nebula 12 + log O/H 
values of 8.51, 8.49, and 8.51 respectively values 
that after adding the fraction of dust tied up in dust 
grains are still smaller than those derived from B 
stars. 

3.3. Comparison of the oxygen abundances of the 
ISM of the solar vicinity with those of the Sun 

and F and G stars of the solar vicinity 

In addition to the evidence presented in section 
2 in favor of large t 2 values, and consequently in fa- 
vor of the O II^l method, there is another indepen- 
dent test that can be used to discriminate between 
the T(4363) method and the O Hrl method that 
consists in the comparison of stellar and H II region 
abundances of the solar vicinity. 

Esteban et al. (2005) determined that the 
gaseous O/H value derived from HII regions of the 
solar vicinity amounts to 12 + log (O/H) = 8.69, 
and including the fraction of O atoms tied up in 
dust grains it is obtained that 12 + log (O/H) = 
8.81 ± 0.04 for the O/H value of the ISM of the so- 
lar vicinity. Alternatively from the protosolar value 
by Asplund et al. (2009), that amounts to 12 + 
log(0/H) = 8.71, and taking into account the in- 
crease of the O/H ratio due to galactic chemical evo- 
lution since the Sun was formed, that according to 
the chemical evolution model of the Galaxy by Ca- 
rigi et al. (2005) amounts to 0.13 dex, we obtain an 
O/H value of 8.84 ± 0.04 dex, in very good agree- 
ment with the value based on the O Url method. 



In this comparison we are assuming that the solar 
abundances are representative of the abundances of 
the solar vicinity ISM when the Sun was formed. 

There are two other determinations of the present 
O/H value in the ISM that can be made from obser- 
vations of F and G stars of the solar vicinity. Accord- 
ing to Allende-Prieto et al. (2004) the Sun appears 
deficient by roughly 0.1 dex in O, Si, Ca, Sc, Ti, 
Y, Ce, Nd, and Eu, compared with its immediate 
neighbors with similar iron abundances, by adding 
this 0.1 dex difference to the solar value by Asplund 
et al. (2009) we obtain a lower limit of 12 + log O/H 
= 8.81 for the local interstellar medium. A similar 
result is obtained from the data by Bensby & Fcltz- 
ing (2006) who obtain for the six most O-rich thin- 
disk F and G dwarfs of the solar vicinity an average 
[O/H] = 0.16; by adopting their value as representa- 
tive of the present day ISM of the solar vicinity we 
find 12 + log O/H = 8.87. Both results are in very 
good agreement with the O/H value derived from 
the Ollflx method. 



3.4. Comparison of the heavy element and helium 
abundances of Ml 7 with those of K dwarf stars 
and with models of galactic chemical evolution 

The best Galactic H II region to determine the 
He/H ratio is M17 because it contains a very small 
fraction of neutral helium and the error introduced 
by correcting for its presence is very small. Ca- 
rigi & Peimbcrt (2008) obtained for M17 a value 
of AY/AZ = 1.97 ± 0.41 for t 2 = 0.036 ± 0.013, 
where Y and Z are the helium and heavy elements 
by unit mass and t 2 was determined observationally. 
By correcting this value considering that the fraction 
of O trapped in dust amounts to 0.12 dex instead 
of 0.08 dex (Mesa-Delgado et al. 2009) we obtain 
AY/AZ = 1.77 ± 0.37. This AY/AZ value is in 
very good agreement with three independent deter- 
minations: the one derived from the chemical evo- 
lution of the Galaxy for the galactocentric distance 
of M17 that amounts to AY/AZ = 1.70 ± 0.4, 
and two AY/AZ determinations derived from K 
dwarf stars of the solar vicinity that amount to 
2.1 ± 0.4 (Jimenez et al. 2003) and to 2.1 ± 0.9 
(Casagrande et al. 2007). On the other hand the 
value AY/AZ = 3.60 ± 0.68 derived from colli- 
sionally excited lines of M17 under the assumption 
of t 2 = 0.00 is not in agreement with the chemical 
evolution models by Carigi & Peimbert (2008) nor 
with the values derived from K dwarf stars of the 
solar neighborhood. 
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TABLE 3 

PRIMORDIAL HELIUM ABUNDANCE 



Method 




Y P 


source a 


HII, {t 2 ± 0.000) 




0.2477±0.0029 


1,2 


HII, (t 2 = 0.000) 




0.2523±0.0027 


1,2 


WMAP + SBBN, (t„ 


= 885.7 ±0.8 s) 


0.2487±0.0006 


3,4 


WMAP + SBBN, (t„ 


= 881.9 ± 1.6 s) 


0.2479±0.0007 


3,4,5,6,7 


WMAP + SBBN, (t„ 


= 878.5 ±0.8 s) 


0.2470±0.0006 


3,5,7 



a (l) Peimbert et al. (2007a); (2) Porter et al. (2005, 2007, 2009); (3) Dunkley et al. 
(2009); (4) Arzumanov et al. (2000); (5) Serebrov et al. (2005, 2008); (6) Mathews 
et al. (2005); (7) Peimbert (2008). 



3.5. Comparison of the primordial helium 
abundance derived from H II regions together 
with that derived from WMAP and SBBN 

In Table 3 we present the determination of the 
primordial helium abundance, Y p , based on two dif- 
ferent methods: a) that based on determination of 
He/H ratios of metal poor H II regions and their ex- 
trapolation to the He/H value for the case of no 
heavy elements, and b) that based on the barion to 
photon ratio derived from WMAP observations to- 
gether with the assumption of Standard Big Bang 
Nucleosynthesis, SBBN. Note that these two values 
might be different if SBBN docs not apply to the 
primordial nucleosynthesis stage of the universe (e.g. 
different number of neutrino families, varying grav- 
itational constant, etc.). 

In Table 3 we present two values derived from H II 
regions (Peimbert, Luridiana & Peimbert 2007a): 
one for constant temperature inside the H II regions, 
that based on the T(4363) value (t 2 = 0.00), and 
the other based on the assumption that the temper- 
ature varies over the observed volume (t 2 ^ 0.00). 
The adopted atomic physics for the helium recombi- 
nation lines is that presented by Porter et al. (2005, 
2007, 2009). We also present three determinations 
based on the WMAP observations by Dunkley et al. 
(2009), and three neutron lifetimes, r„, those by: a) 
Arzumanov et al. 2000, that amounts to 885.7±0.8 s, 
b) Mathews et al. (2005) that amounts to 881. 9± 1.6 
s, and is based mainly on the average of the results by 
Arzumanov et al. (2000) and Serebrov et al. (2005), 
and c) Serebrov et al. (2005, 2008) that amounts to 
878.5 ± 0.8 s. From Table 3 it is clear that the main 
source of error in the WMAP + SBBN method to 
determine Y p is due to the neutron lifetime. 

Moreover it follows from Table 3 that the Y p value 
for (t 2 ^ 0.00) is closer to the Y p value based on 
WMAP and SBBN than the Y p value for (t 2 = 0.00). 
Additional discussion on Y p is presented by Peimbert 



(2008) and by Peimbert et al. (2010). 

4. THE IONIZATION CORRECTION FACTOR 

To derive the total gaseous abundances of a 
given element it is necessary to observe the relative 
amounts of all the ionization stages of that clement 
present in a given nebula. Very often this is not 
possible and one has to correct for the presence of 
the unobserved ions. This correction can be done 
by three different methods: a) an empirical ioniza- 
tion correction curve, b) the use of equations, where 
the correction for the unseen ions is estimated by an 
ionization correction factor based on ratios of other 
observed ions, and c) from tailor made photoioniza- 
tion models. 

Bowen and Wyse (1939) were the first to propose 
the use of what has been called the empirical ion- 
ization distribution curve (e.g. Seaton 1960, Allcr, 
1961, Aller and Lillcr 1968). While PC together with 
other authors started to use equations with ioniza- 
tion correction factors (e. g. Seaton 1968, Rubin 
1969, Peimbert and Torres-Peimbert 1971). 

The calibration of the empirical ionization cor- 
rection factors, ICF's, has been obtained by using 
photoionization models and by observing a given ob- 
ject at different lines of sight with different degrees 
of ionization under the assumption of chemical ho- 
mogeneity. 

The ionization correction factors of N and Ne pre- 
sented by PC have been used widely and are given 
by: 



"M-ICFQl)"^ 



N(0) N(N+) 



N(K) 
and 

AT(Ne) 
A(HT 



A^(H+) N(0+)N(U+Y 



(3) 
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where N(0) = N(0+) + N(0++). 

These two equations present very good approxi- 
mations to the total abundance ratios for those cases 
where the ICF is small, for large values of the ICF 
they should be taken with caution. 

For example in Figure 3 of Peimbert, Torres- 
Peimbert and Luridiana (1995) it can be seen that 
Ne ++ /0 ++ in planetary nebulae increases with de- 
creasing density, indicating that for objects with 
N e < 1000 cm -3 equation (8) is a poor approxima- 
tion to the Ne/O value; this result probably is due 
to the presence of the charge exchange reaction 

0+ 2 + H° - 0+ + H+ 

that permits the coexistence of Nc ++ with + (e.g. 
Hawley & Miller 1977, 1978; Hawley 1978; Pequig- 
not, Aldrovandi & Stasihska 1978; Butler, Ben- 
der & Dalgarno 1979; Pequignot 1980). Ionization 
structure models predict that the lower the density 
the higher the H°/H + ratio, in agreement with the 
charge exchange suggestion and Figure 3. 

Moreover for PNe and H II regions of low de- 
gree of ionization, those with a substantial fraction 
of once ionized O, equation 4 only presents a lower 
limit to the total Ne/H abundance (e. g. Figure 4 
of Torres-Peimbert and Peimbert, 1977; and Figure 
9 of Peimbert, Torres-Peimbert, & Ruiz, 1992). 

Sets of useful ionization correction factor equa- 
tions for gaseous nebulae have been presented by 
many authors (e. g. Kingsburgh and Barlow 1994, 
Izotov et al. 2006). 
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